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Abstract. PKS 1257-326 is a quasar showing extremely unusual, rapid inter- 
stellar scintillation (ISS), which has persisted for at least a decade. Simultaneous 
~ Q ■ observations with the VLA and ATCA, combined with ATCA monitoring over 

I 1 several years, have revealed some properties of the turbulent ionized medium 

responsible for the ISS of PKS 1257-326. The scattering occurs in an unusu- 
ally nearby (~ 10 pc), localized "screen". The scintillation pattern is highly 
anisotropic with axial ratio more than 10 : 1 elongated in a northwest direction 
on the sky. Recent findings and implications for small-scale ionized structures 
in the ISM are discussed. 



1. Introduction 

Interstellar scintillation (ISS) of compact radio sources probes structure in the 
ionized instellar medium of our Galaxy down to extremely small scales (~ 
10 5 km). PKS 1257-326 (Galactic coordinates I = 305.2°, b = 29.9°) is a flat 
spectrum, radio-loud quasar at z = 1.256 which ex hibits intrahour flux density 
variability at frequencies of several GHz due to ISS ( Bignall et al.ll2003l ). It was 



the third confirmed "intrahour variable" (IHV) quasar, after PKS 0405—385 
dKedziora-Chudczer et~al"1ll997T ) and J1819+3845 (jDennett-Thorpe fc de Bruvn 



20001 ). 

The recent large-scale Micro Arcsecond Scintillation- Induced Variability (MA- 
SIV) VLA Survey at 5 GHz showed that while a large fraction of flat-spectrum 
radio sources exhibit variations of typically a few percent on timesc ales of or- 



der a day, the IHV phenonemon is extremely rare (jLovell et al.ll2003l . and these 



proceedings). Only a handful of quasars are known to show large variations on 
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timescales of a few hours or less. While a source must be compact in order to 
scintillate, short scintillation time-scales have been attributed to nearby scatter- 
ing screens within a few te ns of parsec from the Sun (jDennett-Thorpe fc de Bruvn 
20001 : iRickett et al.ll2002h . 



2. Observations and analysis 

We began monitoring PKS 1257-326 with the ATCA at 4.8 and 8.6 GHz in 
2000 following the discovery of rapid variability ( Bignall et al.l [20031 ). We found 



an annual cycle in the timescale of variability which is evidently consistent over 
several years of observations at both frequencies. The annual cycle results from 
the change in scintillation velocity due to the Earth's orbital motion. Further- 
more, simultaneous observations of PKS 1257—326 with the VLA and ATCA 
revealed large delays of up to 8 minutes bet ween the variability pattern arrival 
times at each telescope ( Bignall et al. 20061 ) . as shown in Figure [TJ The two- 



station time delays were measured in three separate epochs, simultaneously at 
two frequencies and on two consecutive days in eac h epoch. Unlike for the case of 



J181 9+3845 observed with the VLA and WSRT (lDennett-Thorpe & de Bruvn 



2002), where the projected baseline rotated through a large angle during the 



simultaneous observation and the time delay changed sign, for PKS 1257—326 
observed with the VLA and ATCA we detect no significant change in delay over 
the course of an observation. The mean pattern arrival time delays for each 
epoch were 483 ± 15 s in 2002 May, 333 ± 12 s in 2003 January, and 318 ± 10 s 
in 2003 March. The observed annual cycle and time delays not only demon- 
strate conclusively that the rapid, large-amplitude variability of PKS 1257—326 
is entirely due to ISS, but moreover these results can be used quantitavely to 
determine the bulk velocity of the scattering plasma as well as the characteristic 
length scale of the scintillation pattern in two dimensions. 

In our analysis, the scintillation pattern is assumed to be elliptical with 
minor axis characteristic scale length ttmin? 

and major axis a ma j = i?a m j n oriented 
along the vector S = (cos /3, sin (3). The scintillation velocity is written in the 
form v(T) = viss — Vm (T) _ (v n (T), v a (T)), and v aries on an annual cycle as 



a function of time, T ( Macquart Jauncevi I2OO2I ) . Here we define Vq(T) as 



the velocity of the Earth with respect to the solar system barycenter. In this 
notation, the ti me delay expected betwe en two telescopes displaced by a distance 
r = (r a ,r 5 ) is (|Coles fc Kaufmanlll978T ) 



At = r-v+(fi 2 -l)(rxS)(vxS) 
v 2 + {R 2 - l)(v x S) 2 



where the dependence of v on T is suppressed. The scintillation timescale is 

, 0" m i n 



v 2 + {R 2 -l)(v x S) 2 



(2) 



We define i S cint as the point where the intensity autocorrelation function decays 
to 1/e of its maximum value. 

In practice neither the annual cycle nor the time delay experiment datasets 
suffice to uniquely determine all five scintillation parameters, viss a, viss <5> ft, R 
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PKS 1257-326, 2002 May 14, 4.9 GHz 
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Figure 1. Simultaneous flux density measurements of PKS 1257—326 at the 
VLA and ATCA on 2002 May 14 at 4.9 GHz. The pattern takes 8 minutes 
to drift over the 10 4 km baseline, yet the annual cycle in timescale implies a 
scintillation velocity significantly larger than 20 km s" 1 at this time of year. 
The large observed time delays in May are due to a highly anisotropic pattern 
elongated almost perpendicular to the baseline and moving with a significant 
velocity component parallel to the baseline. 



and a m i n , at any one observing frequency, so we fit to both the time delay and the 
annual cycle data simultaneously. We also fit both frequencies simultaneously, 
imposing the additional constraint that the scintillation velocity is identical at 
the two frequencies; R and (3 at the two frequencies are still assumed indepen- 
dent. 

2.1. Anisotropy 

When the scintillation pattern is highly anisotropic, there are degenerate solu- 
tions for the axial ratio R and the component of velocity parallel to the major 
axis. However the characteristic length scale along the minor axis, a m i n , can still 
be uniquely determined. In weak scattering, as observed for PKS 1257—326 at 
frequencies of ~ 5 GHz and above, the minor axis length scale is related to the 
Fresnel scale, rp = v 'cL j '(2nv), where v is the observing frequency and L is the 
scattering screen distance. Thus a m i n can be used to constrain L. 

For PKS 1257—326 we find that a highly anisotropic scintillation pattern is 
required to fit the data. The position angle of anisotropy is well constrained to 
within a few degrees but there is a degeneracy between the two components of 
the scattering screen velocity, and the pattern axial ratio R is poorly constrained. 
Figure [2] shows the timescale data and annual cycles from several models impos- 
ing different constraints on the parameters. For the fits shown, we weighted x 2 
to give the (comparatively few) time delay measurements equal importance in 
the fit relative to the characteristic timescale measurements. The dotted line in 
Figure [5] show that assuming an isotropic scintillation pattern results in a very 
poor fit to the data. 

Based on the similarity of the patterns observed at the two telescopes, 
we can determine a lower limit for R. When the scintillation velocity has a 
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PKS 1257-326, ATCA 
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Figure 2. Annual cycle in characteristic timescale, i S cint, observed in ATCA 
data at 4.8 GHz between 2000 and 2003. The plotted lines show model annual 
cycles for best fits to the combined time delay and annual cycle data with (a) 
the scintillation pattern fixed to be isotropic (R = 1); (b) no constraints on 
the parameters; (c) pattern axial ratio limited to R < 12; (d) screen velocity 
fixed to the local standard of rest. 



significant component perpendicular to the baseline r, as is the case for the 
May time delay data, we can expect a small amount of spatial decorrelation of 
th e scintillation patter n over the ATCA-VLA projected baseline. As detailed 
in iBignall et al.l ( 20061 ) after careful processing of the data we can constrain 



the axial ratio R to be at least ~ 12 : 1 based on the degree of correlation 
between the patterns seen at each telescope. Figure [3] shows the annual cycle in 
scintillation velocity for a scattering screen moving at — 49. 2 km s -1 in RA and 
11.5 kms -1 in Dec with respect to the Sun. The position angle of anisotropy 
in the scintillation pattern (—55° North through East) is also illustrated by the 

filled ellipse in Figure [31 

As discussed by iRickett &: Colesl ( 20041 . and these proceedings), there are 



several lines of evidence for anisotropic turbulence in the ISM. Moreover the 
intensity autocorrelation functions for PKS 1257—326 disp lay strong oscillator y 
behaviour which is a signature of anisotropic scattering ( Rickett et al.l [20021 ) . 



Therefore it seems likely that the large anisotropy in the scintillation pattern re- 
sults from anisotropic scattering rather than anisotropic source structure. Highly 
anisotropic scattering suggests a plasma controlled by a magnetic field and lo- 
calized in a relatively thin "screen" . 

2.2. Distance to the scattering screen 

Our best fit models give a m i n = 4.24 ± 0.08 x 10 km at 5 GHz, and a m ; n = 
3.43 ± 0.06 x 10 4 km at 8.6 GHz. For anisotropic scattering, our analysis shows 
that our definition of a m i n corresponds to 0.78rp for a point source in weak 
scattering. If the source angular size is larger than the angular Fresnel scale 



Scintillating Quasar PKS 1257-326 5 




80 60 40 20 

V a (km s- 1 ) 



Figure 3. Scintillation velocity projected onto the plane of the sky. The 
large ellipse shows the annual variation of the Earth's velocity with respect to 
a scattering screen moving at —49.2 km s^ 1 in RA and 11.5 km s _1 in Dec with 
respect to the Sun, corresponding to solution (c) in Figure O Dotted lines 
indicate the velocities on the days of the time delay observations. The 1-cr 
contour of the best fit solution for fiss.a versus viss,s is plotted around the best 
fit screen velocity at the center of the large ellipse. The component of velocity 
parallel to the long axis of the scintillation pattern is poorly constrained. If 
the scattering screen were moving with the local standard of rest, the large 
ellipse would be centered at position x marked in the top right-hand corner. 
The black filled ellipse represents the anisotropy (position angle and axial 
ratio) of the best fit scintillation pattern at 4.9 GHz. 



then a m j n will be larger than rp . Thus our data put an upper limit on the screen 
distance of L <^ 10 pc. 

2.3. The influence of source structure? 

Although there is evidence for the scattering medium being largely responsible 
for the high degree of anisotropy in the scintillation pattern of PKS 1257—326, 
it is possible that source structure also plays a role in the observed scintillation. 

If the scintillating component is circularly symmetric, then assuming a 
screen distance of 10 pc, approximately 100 mJy would have to be contained 
within an angular size of ~ 30 /ias in order to display the observed large modu- 
lations. This implies a brightness temperature of at least ~ 10 13 K, requiring 
a Doppler factor of 10-20 to reduce the intrinsic brightness temperature below 
the inverse Compton limit. Alternatively, if the "core" happened to be elon- 
gated, or consist of several components aligned in the same direction as the ISM 
anisotropy, then this would reduce the brightness temperature requirements. 
While an alignment of source structure with ISM anisotropy seems an unlikely 
coincidence, it is conceivable that t his could be another "selection effect" for 
IHV. Both VLA (Big nall et aT1l2006ft and VLB A (lAkaba. Bignall Revnoldsl 



2006, in press) images of the source show broad "jet" components extended to 



the northwest of the core. Speculatively, if the unresolved VLBI core consists 
of several components on sub-mas scales also aligned in this direction, then this 
implies an approximate alignment of the jet with the major axis of anisotropy 
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in the scintillation pattern. In this case, since there is effectively only significant 
scattering in the direction perpendicular to the jet, the scattering screen would 
not "resolve" structure along the jet axis. 

3. Conclusions 

Simultaneous observations with the VLA and ATCA, combined with ATCA 
monitoring over several years, have revealed some properties of the scattering 
plasma responsible for the rapid, intrahour-timescale scintillation of PKS 1257- 
326. The scattering occurs in an unusually nearby (~ 10 pc), localized "screen". 
Based on the longevity of the scintillation we can derive a lower limit of 10 15 cm 
on the linear extent of the screen. Intrahour variability is extremely rare, sug- 
gesting that these nearby screens cover a tiny fraction of the sky. The tur- 
bulence is highly anisotropic, producing a scintillation pattern with axial ratio 
more than 10 : 1, elongated in a northwest direction on the sky. It was suggested 
by J. Linksy (private communication, this conference) that the scattering could 
possibly be produced in turbulent regions of interaction at the edges of local 
clouds. It will be important to investigate the possible ISM structures respon- 
sible for the rapid scintillation of some quasars, not only to better understand 
the tiny-scale structure in the ionized ISM, but also to gauge the effects of ISS 
on observations with future sensitive cm-wavelength telescopes. 
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